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to Abstract 

An experiment was conducted at West Lafayette, Indiana in 1979 to quanti- 
tatively describe the Interaction of the solar illumination angle and row 
azimuth angle on the measured reflectance factor (RF) of soybean canopies 
consisting of IX plots. Nine of the plots were planted in 71 cm wide rows; 
the other plots were of bare soil and soybeans with 100 percent soil cover. 
Reflectance factor data in four spectral bands, 0. 5-0.6, 0.6-0. 7, 0.7-0. 8, and 
0. 8-1.1 ym, were taken at 15 minute intervals during three clear days, August 12 
and 31 and September 19 over nine plots of differing azimuthal direction with a 
Landsat-band radiometer (Exotcch Model 100) at 5.2 meters above the soil. 

Diurnal changes of nearly 140 percent were observed in the red wavelength 
region when canopies covered 64 percent of the soil. The amount of shadow 
observed was a function of the plant geometry and row width. As soil cover 
approached 100 percent, the diurnal changes diminished. A function that 
describes the solar illumination angle with respect to the row azimuth 
explained most of the diurnal variation in the measured RF. Variation in near 
infrared response was much less and did not appear to be as strongly related 
to sun-row angle interactions. The ratio, near Infrared/red, was highly sen- 
sitive to sun angle-row direction interactions, whereas the greenness function, 
utilizing all four spectral bands, was not. 
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Introduotion 

Understanding the effect of the interactions between solar 
illumination and crop canopy geometry on the spectral response is 
necessary to utilize reflectance factor data effectively. Numerous 
models have been proposed bo explain and predict the measured 
reflectance factor of plant canopies as a function of plant geometry, 
sun angle, and view angle (Suits, 197(?| Smith eb al., 1975; Richardson 
eb al., 1975). The models by Suits and Smith deal with a canopy with no 
horizontal spatial variations. 

Richardson eb al. (1975) modeled the reflectance of a row crop, 
with distinct horizontal spatial variations, as a function of plant, 
soil, and shadow components. By illuminating a surface covered with 
various shaped objects, Egbert (1977) was able bo explain 80 bo 85 
percent of the variance in the reflectance measurements due to shadows. 
A .lodel suggested by Jackson et al. (1979) assumes on incomplete canopy 
of rectangular-shaped rows. The fractions of sunlit and shaded soil and 
vegetation viewed are calculated as a function of view angle for a 
particular canopy condition, described by plant cover, helght/wldbh 
ratio, row spacing and direction, time of day, day of year, latitude, 
and size of the radiometer resolution element. 

Studies of the effect of sun zenith angle on reflectance generally 
have supported the predictions o. the Suit's canopy reflectance model 
that the refleobanoo factor should increase as the solar elevation 
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increases (Colwell, Chance and LeMaater, 1977). Colwell (197*0 
attributes this to changes in the amount of shadow within the canopy. 
Field data have shown minor to significant Increases in the infrared 
response with decreasing sun elevations (Puggin, 1977; Chance and 
LcMaster, 1977; Jackson et al., 1979). Crecelius (1978) noted symmetric 
and non-symmetric components about solar noon that influenced the 
observed variation in reflectance throughout the day. The symmetric 
component, solar angle, explained the majority of the observed 
variation. Other effects, such as drying of the soil surface and plant 
wilting will be asymmetric about .solar noon and may be significant 
factors to consider. 

Further investigation of reflectance factor data taken in 1978 over 
incomplete soybean canopies revealed possible time of day effects in the 
Landsat band regions as illustrated in the red, 0.6-0. 7 pm* and the near 
infrared, 0, 8-1.1 pm, in Figure 1. Plots were planted in a north-south 
row direction. Both bands were plotted with and without a 1.5 hour time 
restrictic-n about solar noon. Low responses were noted over those plots 
that were measured more than 1.5 hours from solar noon. These changes 
in RF resulted from shadows between the rows and a lower response from 
the soil component. 

The objective of this research was to model the reflectance of a 
soybean canopy as a function of row width, row direction, and solar 
azimuth and zenith angles. The hypothesis was that by varying only the 
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row direction, the variation in reflectance would be explained entirely 
by changes in sun Mnibh and azimuth angle with respect to row 
direction. 


Materials and Methods 


Experimenta l Co nditions 

The experiment was conducted in 1979 on the Purdue Agronomy Farm. 
Soybeans ( Glycine max (L.) Merr. "Amsoy 71") were planted on a Chalmers 
silty clay loam (typic Argiaquoll) on June 25, 1979. Soybeans ware used 
because they have dense foliage with distinct row patterns through much 
of the season. This is in contrast bo many of the other major crops 
such as wheat ard corn that have a much more complex canopy geometry and 
shadow pattern. Because of extended periods of cloudy days early in the 
season, spectral data were collected only on development stages after 
full bloom. 

The experiment consisted of 11 randomly arranged plots which were 
3.5 m wide and 5.P m long (Figure 3). Nine plots were planted in 71-om 
rows with the following azimuthal directions: 0-180, 30-210, 60-2^0, 
90-270, 105-285, 120-300, 135-315, 150-330, and 165-345 degrees from 
north. Another plot was planted in 26-cm wide east-west and north-south 
rows to obtain a canopy with negligible row effects. A bare soil plot 
was included to monitor the sunlit soil background reflectance of the 
soybean plots. Row directions were selected to favor data collection 
during the morning hours when cloud-free conditions were more likely. 
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Three development stages with 65 » 78, and 9^ percent soil cover on 
the 71-cm rows were represented with the three measurement dates. The 
canopy with 78 percent soil cover was obtained by trimming a near full 
canopy just prior to the start of senescence. The cross sectional shape 
of the canopy was determined by placing a large piece of poster board in 
the canopy, perpendicular to the row azimuth, at several locations and 
drawing the perimeter of the canopy on the board. The canopy shapes for 
each date are illustrated in Figure 3. 

Spectral Measurements 

Radiance measurements, used bo determine reflectance factor (RF), 
were taken over all the plots with a Landftat band radiometer (Exotech 
Model 100) at 15-minute intervals throughout the day on three clear days 
(August 12, August 31? and September 19). Nicodemus et al. (1977) and 
Robinson and Biehl (1979) desoibe the conditions and proceedures for 
obtaining the reflectance factor, which closely approximates the 
bidirectional reflectance factor. The Exobech 100 is a ^1-band 
radiometer with a 15-degree field of view that acquires data in the 
following wavelength regions: 0.5-0. 6, 0.6-0. 7, 0.7-0. 8, and 
0.8- 1.1 pm. Data were taken only under near cloud-free conditions 
(especially in the vicinity of the sun). 

A mobile truck-mounted radiometer system was used for quick and 
efficient data collection in the field. A boom mounted on the back of 
the truck permitted the radiometer and a motor-driven camera to be 
placed 5.2 m above the crop canopy and 3.5 m from the truck. Spectral 
data were collected over two locations in each plot on August 12 and 


Height (cm) 
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over four locations In each plot on August nnd September 19 • The 
instruments were carefully leveled to obtain all spectral data at a 
nadir look angle* Several measurements were taken over each plot and 
averaged to insure a representative sampling of the plot and bo avoid 
biased values for on-row or off-row measurements. Measurements in all 
bands were taken concurrently and recorded by a printing data logger. 
During data oolleobionf photographs were taken periodically over each 
plot for soil cover determination and shadow assessment. 

Agronomic Measurements 

Agronomic measurements included plant height i leaf area index, 
maturity stage (Fehr and Caviness, 1977), surface soil moisture, total 
fresh and dry biomass, and stem, pod, and green leaf dry biomass. 
Percent soil cover was determined by placing a grid over the verbicfil 
photograph and counting the intersections occupied by green vegetation. 

Data Analysis 

The reflectance factor data were analyzed as band means. The 
reflectance data were transformed into greenness as described by Kauth 
and Thomas (1976) for Landsat MSS data and modified for speotroraoter 
data (Malila and Gleason, 1977). The data transformation was; Greenness 
= [(BandB » 0.17289) + (Band.Jl » 0.59538)] - [(Bandl » 0.i|8935) + (Band2 
» 0.6l2i19)j . Bandl to band4 refer to the RF measured in the four 
Landsat bands. The near infrared/red reflectance ratio |^(0. 8-1.1 
um)/(0,6-0*7 pm)j was also considered in the analysis. Analysis of 
variance and Newman-Keuls tests were performed to determine significant 
effects of row-solar angle interaction and RF. 
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Results find Dlsousslon 

The maximum response of RF to changes in sun angle occurred when 
the sun azimuth angle was equfti to the row azimuth angle (Figure Jla). 
Diurnal changes in RF of nearly 1*10 percent were observed in the red 
wavelength region, 0.6-0. 7 )jm, on August IF. The highest reflectance 
values were obtained when the soil was sunlit and the lowest, when the 
soil was shaded. Diurnal variations in the RF .in the near Infrared 
wavelength region, 0. 8-1.1 pm, were lower (relative to the minimum RF 
value observed) than that noted in the visible region and not as clearly 
related to sun-row interactions (Figure *}b). Note the absolute changes 
in RF are about the same. The shadows of the near infrared region may 
nob be as dark as those observed in the visible region due to low 
pigment absorption and multiple scattering in the canopy (Colwell, 
197*0. 

The effect of sun-row azimuth interactions are shown in Figure 5. 
The reflectance was plotted over time for three plots of different row 
directions. The peak response in the red wavelength region for the 
three plots was not only at different times, but also in order of the 
row azimuth. Again the peak response was when the sun was shining down 
the rows, lighting the soil surface, and thus giving a higher 
reflectance reading. 

t 

The diurnal response in the red wavelength region for two of the 
key canopy components, sunlit soil and vegetation, are shown in Figure 
6. Very little change in reflectance factor was observed as a function 
of zenith angle for either the plot containing bare soil or the plot 
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Solar Azimuth - Row Azimuth 


Figure 4. Changes in Lhe BRP in (a) the red wavelength band (0.6-0. 7 
;im) and (b) the near infrared wavelength band (0. 8-1.1 jjm) plotted 
against the difference between solar and row azimuth. Row azimuth*' 1 80°. 
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Figure 5. RF in the red wavelength region (0.6-0.7 pm) for three 
row directions over time on August 12, 1979. 
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Figure 6. Reflecconcc factor in the rod v«ivolongt!» (0.6-0. 7 ^im) Cor the 
bare v^oil and full vegetative ennopy plots against changes in zenith angle 
on August 12, 1979. 
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with 100 percent solX cover. Note the Inrge Ulfferences between the 
sunlit soil and sunlit vegetation. Interaction of the canopy shape and 
size, and soil width with the aim angle produces varying amounts of 
shadow east on both the soil and the vegetation. Thus, it may be that 
the diurnal variations in RF observed in Figures and b were caused by 
changes in the amount of shadow in the field of view of the instrument. 

Equations to predict the shadow cast by rectangular or spherical 
rows as the sun zenith and azimuth change throughout the day have been 
defined by many investigators (Idso and Baker, 197R| Jackson et al., 
1979; Verhoff and Bunnik, 1978). For this study, an equation was used 
bo express the solar zenith (0) and azimuth angle (i|t) in relation to a 
projected ray onto a plane perpendiouXar to the row azimuth. This 
function, called the projected solar angle, 9sp = tan~1 (ban9aln(|i) , is 
illustrated in Figure 7* 

The response of the red, near infrared, near infrared/red 
reflectance ratio, and greenness bransformation bo changes in 0sp on 
August 12 are presented in Figure 8. The near infrared/red ratio and 
the greenness transformation are often related bo changes in plant 
biomass, soil cover, and/or leaf area index. If the diurnal changes are 
expressed as a percent increase in response relative to the minimum 
value observed that day, the red and near infrared/red ratio were quite 
sensitive bo changes in 0sp, whereas the greenness and the near Infrared 
region were nob. The near infrared wavelength region does show similar 
absolute changes in response during the day, but the pattern is nob as 
clearly related to the changes in sun-row angle and the variation about 
the mean for any given Ssp is muoh higher for the near infrared region. 



Date 

Time 

Solar Angle 
Zenith Azimuth 

®sp 

Aug 15 

12100 

26,3 

180,0 

26.3 

Aug 21 

2:09 

40.1 

234.2 

26,3 

Sept 10 

4:49 

73.7 

261,7 

26,3 


Figure 7. Xlluscraclcm of the solar projected angle 0 as observed 
when looking to the west. Tliree dates and the correspon§Sng observation 
times result In the same shadow pattern and 6 
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The response in the red (0.6-0. 7 pm) wavelength band has been 
plotted as a function of 0sp for the three diurnal studies in 1979 
(Figure 9). The canopies with lower soil covers, 64 and 78 percent, 
showed greater changes in reflectance due to changing sun angle than the 
near full canopy of 94 percent soil cover. The RF of the canopies with 
94 percent soil cover changed only slightly more during the day than the 
RF of the full canopy. The first two dates appeared to have two 
functions present, the first being highly dependent on <j)sp, the second 
independent of 0sp. 

The dependent zone, where the RF is changing rapidly with changes 
in 0sp, is a function of the sunlit soil reflectance and the vegetation 
reflectance (Figure 9). The variation about the mean might be due to 
local variations in soil cover or, possibly, instrument position about 
the row at low sensor altitudes. Some of this variation was thought bo 
be due bo the interaction of sun zenith angle with the surface roughness 
of the canopy, with large zenith anglen causing longer shadows and thus 
lower reflectance. However, no evidence of this was apparent from the 
analysis of the data. 

In the independent region (Figure 9), where the soil surface was 
completely shadowed, the measured reflectance was a function of by one 
variable, the percent soil cover. Just as for the dependent zone, local 
variations in soil cover might cause the observed amount of variation 
about the mean. The critical angle, beyond which a change in the 
projected angle no longer results in a change in RF, shifts to lower 
0sp's for higher soil covers or canopy heights. 
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Figure 9. Relationship between RF in the red wavelength band (0.6- 
0.7 1 m) and projected solar angle Osp for three different canopies. 
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Sumniary and Conoluaions 

The objective of this research was to identify the physical causes 
of diurnal changes in the reflectance factor of a row crop canopy. This 
is important when the measured reflectance factor in the visible region 
of a given plot may vary 100 percent or more during the day due to 
varying amounts of shvidow within the canopy. A function called the 
projected solar angle, Osp, that includes both the solar zenith and 
azimuth angle plus the row azimuth angle, descibes the changes in shadow 
and, thus, the diurnal changes observed in reflectance factor. Canopy 
geometry was a key factor determining both the diurnal range of the 
spectral response and the critical angle where further increases in Osp 
did not lower the RF measured. The soil component of the soene was then 
completely shadowed. 

The effect of solar zenith angles between 20 and 60 degrees on the 
measured reflectance was found to be nonsignificant when the RF was 
measured at nadir over all the plots, including the bare soil and full 
canopy plots. The near infrared RF and greenness function were not as 
sensitive to changes in solar illumination angle in the row orop canopy 
observed as the visible region and near infrared/red reflectance ratio. 
These variables may thus prove to be useful in relating spectral 
response to such agronomic variables as percent soil cover, leaf area 
index, and plant biomass over a wide range of illumination angles. 

The results indicate that changes in canopy shadowing may be a 
significant factor, particularly in the visible wavelength region, 
influencing the spectral reflectance of crop canopies. A physical model 
accounting for this variation was developed. It will be used in future 
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investigations to simulate the variation which may be expected with 
varying row direction, amounts of canopy cover, date, time of day, and 
latitude. 

Future studies should include a wider range of solar zenith and 
azimuth angles and more row azimuth angles. This objective, along with 
a decrease in plot to plot variability, could be obtained by placing the 
plot to be studied on a turntable. This would allow for a quick change 
in the row direction with a minimum of plot to plot variation. To study 
effects due to the solar zenith angle on soybean row crops more 
effectively, measurements should be taken at low latitudes, where the 
range in zenith angles will be the greatest. 
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